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Abstract 

We present a comparison of electron-phonon interaction in NbBi and MgBi-, cal- 
culated using full-potential, density- functional-based methods in P6/mmm crystal 
structure. Our results, described in terms of (i) electronic structure, (ii) phonon 
density of states F{ui), (iii) Eliashberg function a 2 F(uj), and (iv) the solutions of 
the isotropic Eliashberg gap equation, clearly show significant differences in the 
electron-phonon interaction in NbBi and MgBi ■ We find that the average electron- 
phonon coupling constant A is equal to 0.59 for MgBi and 0.43 for NbBi, leading 
to superconducting transition temperature T c of around 22 K for MgBi and 3 K 
for NbB 2 . 



The lack of success in finding superconductivity in other diborides with super- 
conducting transition temperature, T c , close to that of MgBi [1] underscores 
the complex nature of interaction responsible for superconductivity in MgB 2 . 
In MgB 2 the complexity is further compounded by the presence of multi- 
faceted Fermi surface [2,3] and a highly anisotropic electron-phonon coupling, 
A(k, k'), over the Fermi surface [4,5]. The dependence of superconducting prop- 
erties on such details has ensured that we do not know, as yet, the exact nature 
of interaction leading to superconductivity in MgB 2 . 

Within Eliashberg-Migdal theory [6,7] of superconductivity, a reliable descrip- 
tion of the superconducting state requires an accurate knowledge of A(k, k') 
and the renormalized electron-electron interaction, /x*, which are used as input 
to the fully anisotropic gap equation. The present computational capability 
allows us to evaluate A(k, k') accurately using density-functional-based meth- 
ods but, unfortunately, /x* cannot be evaluated. However, it is reasonable to 
assume that /x* varies between 0.1 to 0.2 [4,5]. Thus, the electron-phonon cou- 
pling A(k, k'), which is a normal state function, must contain signatures of 
superconducting state. 
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In an attempt to identify some of the unique features of electron-phonon in- 
teraction in MgB 2 vis-a-vis other diborides we have studied (i) the electronic 
structure, (ii) the phonon density of states (DOS), (iii) the Eliashbrg function, 
and (iv) the solutions of the isotropic Eliashberg gap equation for NbB 2 and 
MgB 2 in P6/mmm crystal structure. 

The choice of NbB 2 has been motivated by the recent reports [8,9] of supercon- 
ductivity, with T c going up to 9.2 K, under pressure in hole-doped Nb x B 2 . Ear- 
lier experiments have shown superconductivity in stoichiometric NbB 2 [9,10] 
as well as Boron-enriched NbB 2 [11] samples. The reported T c for stoichiomet- 
ric NbB 2 varies from 0.62 K [10] to 5.2 K [9], while for Boron-enriched NbB 2 .5 
the T c is found to be 6.4 K [11]. We also note that, recently, Kaczorowski et 
al. [12] did not find any superconductivity in NbB 2 down to 2K. 

We have calculated the electronic structure of NbB 2 and MgB 2 in PQ/mmm 
crystal structure with optimized lattice constants a and c, as given in Table I. 
The lattice constants a and c were optimized using the ABINIT program [13] 
based on pseudopotentials and plane waves. For studying the electron-phonon 
interaction we used the full-potential linear response program of Savrasov 
[14,15], and calculated the dynamical matrices and the Hopfield parameter. 
These were then used to calculate the phonon DOS, F(u>), the electron-phonon 
coupling A(k, k'), and the Eliashberg function, a 2 F(u), for NbB 2 and MgB 2 . 
Subsequently, we have numerically solved the isotropic Eliashberg gap equa- 
tion [6,7,16] for a range of ji* to obtain the corresponding T c . 

Based on our calculations, described below, we find significant differences in 
the phonon DOS and the Eliashberg functions of NbB 2 and MgB 2 . In par- 
ticular, we find that the average electron-phonon coupling constant is equal 
to 0.59 for MgB 2 and 0.43 for NbB 2 , leading to superconducting transition 
temperatures of around 22 K for MgB 2 and 3K for NbB 2 . 

Before describing our results in detail, we provide some of the computational 
details of our calculation. The structural relaxation was carried out by the 
molecular dynamics program ABINIT [13] with Broyden-Fletcher-Goldfarb- 
Shanno minimization technique using Troullier-Martins pseudopotential [18] 
for MgB 2 and Hartwigsen-Goedecker-Hutter pseudopotential [19] for NbB 2 , 
512 Monkhorst-Pack [20] Appoints and Teter parameterization for exchange- 
correlation. The kinetic energy cutoff for the plane waves was 110 Ry for MgB 2 
and 140 Ry for NbB 2 . The charge self-consistent full-potential LMTO [14] 
calculations were carried out with the generalized gradient approximation for 
exchange-correlation of Perdew et al [21] and 484 Appoints in the irreducible 
wedge of the Brillouin zone. For MgB 2 , the basis set used consisted of 3k 
panels and s, p, d and / orbitals at the Mg site and s, p and d orbitals at the 
B site. In the case of NbB 2 , we included 2k panels and s, p and d orbitals at 
the Nb site. In all cases the potential and the wave function were expanded 
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Table 1 

The calculated lattice constants a and c. The experimental lattice constants for 
MgB 2 [17] and NbB 2 are shown in the parentheses. 





a (a.u.) 


c (a.u.) 


MgB 2 


5.76 (5.834) 


6.59 (6.657) 


NbB 2 


5.81 (5.837) 


6.10 (6.245) 



Table 2 

The site- and /-resolved electronic densities of states, in st/ (Ry — atom), at the Fermi 
energy in MgB 2 and NbB 2 calculated at the optimized lattice constants using the 
full-potential LMTO method. 



alloy 


element 


s 


P 


d 


/ 


MgB 2 


Mg 


0.47 


0.72 


0.94 


0.08 




B 


0.06 


3.36 


0.15 




NbB 2 


Nb 


0.02 


.09 


9.54 






B 


0.07 


1.34 


0.19 





up to Imax — 6. The muffin-tin radii for Mg, B, and Nb were taken to be 2.4, 
1.66, and 2.3 atomic units, respectively. 

The calculation of dynamical matrices and the Hopfield parameters for MgB 2 
were carried out using a 6 x 6 x 6 grid while for NbB 2 we used a 4 x 4 x 4 
grid resulting in 28 and 12 irreducible q-points, respectively. For Brillouin 
zone integrations in MgB 2 we used a 6 x 6 x 6 grid while for NbB 2 we used 
8x8x8 grid of k-points. The Fermi surface was sampled more accurately 
with a 24 x 24 x 24 grid of k-points using the double grid technique as outlined 
in Ref. [15]. 

Here, we like to point out the reasons for carrying out the linear response 
calculation for MgB 2 in spite of earlier calculations by Kong et al. [4] and 
Choi et al. [3,5]. The linear response calculation by Kong et al. is similar 
to the present approach, while Choi et al. used pseudopotentials and frozen 
phonon method to evaluate the electron-phonon coupling A(k, k'). The present 
approach differs from the work of Kong et al in the selection of q-points 
and the Brillouin zone integrations. As a result Kong et al. find the average 
electron-phonon coupling constant A = 0.87 ± 0.05, which is much higher 
than the value of A = 0.61 as reported by Choi et al. [3,5], as well as the 
experimental values of 0.58 [22] and 0.62 [23] as deduced from specific heat 
measurements. Our calculated value of A is equal to 0.59, in close agreement 
with the work of Choi et al. 

A comparison of site- and /-resolved electronic density of states of NbB 2 
[24] and MgB 2 [25] at the Fermi energy is given in Table II. A further 
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Fig. 1. The phonon density of states F(u) of MgB 2 and NbB 2 calculated using the 
full-potential linear response method as described in the text. 

decomposition of the densities of states (st/Ry) in terms of cubic harmon- 
ics reveals the dominance of B p electrons at the Fermi energy in MgB 2 
(Px( y ) — 0.9, p z = 1.55) than in NbB 2 (p x ( y ) = 0.43, p z = 0.48). In addition, 



in NbB 2 , the Nb <i-electrons (d xy ( x 2_ 



-y 2 ) 



1.83, d yz ( zx ) = 1.66, d 3 2 2 -i = 2.56) 



are present in substantial amount, indicating a more active role for Nb in 
determining the possible superconducting properties of these materials than 
played by Mg in MgB 2 . 



In Fig. 1 we show the phonon DOS F(u) of NbB 2 and MgB 2 calculated using 
the full-potential linear response program as described earlier. For MgB 2 we 
can clearly identify four significant peaks in the phonon DOS at 33, 53, 79 
and 96 meV, respectively. The peak at 33 meV is related to the van Hove 
singularity in the acoustical mode [5], and it involves the motion of Mg atom 
and B atoms separately. However, the region around the peak at 53 meV 
results from the motion of both Mg and B atoms. The phonon DOS around 
79 and 96 meV peaks are due to the coupled motion oiB — B atoms in the x—y 
plane. In particular, the peak at 79 meV corresponds to the in-plane B — B 
bond stretching mode, and in Ref. [5] it is located at 77meV. Similarly, for 
NbB 2 , the peak in F(u) at 32 meV is dominated by the motion of Nb atom, 
while the region around 65 — 70 meV results from the coupled motion of Nb 
and the two B atoms. Not surprisingly, 81 meV peak in the phonon DOS of 
NbB 2 corresponds to the in-plane B — B motion. In contrast with the phonon 
DOS in MgB 2) the phonon DOS in NbB 2 around 106 meV results from the 
displacements of both Nb and the two B atoms. 

To see the strengths with which the different modes of the ionic motion couple 
to the electrons, and thus are capable of influencing the superconducting prop- 
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Fig. 2. The Eliashberg function a 2 F(uj) of MgB 2 and NbB 2 calculated using the 
full-potential linear response method as described in the text. 

erties the most, we show in Fig. 2 the Eliashberg function a 2 F(uu) of NbB 2 
and MgB 2 calculated as described earlier. The most striking feature of Fig. 2 
is the overall strength of the electron-phonon coupling in MgB 2 as compared 
to NbB 2 . We find that the average electron-phonon coupling constant A is 
equal to 0.59 for MgB 2 and 0.43 for NbB 2 , which clearly shows that MgB 2 is 
more likely to show superconductivity with a higher T c than NbB 2 . 



Further analysis of the Eliashberg function, as shown in Fig. 2, reveals the 
importance of the in-plane B — B bond-stretching optical phonon mode in 
MgB 2 , which gives rise to the dominant peak at 79meV. The other peaks in 
the phonon DOS of MgB 2 , such as the peaks at 33 and 53 meV, couple weakly 
with the electrons at the Fermi energy. Thus the motion of Mg atom plays 
a relatively insignificant role in determining the superconducting properties 
of MgB 2 . In contrast, in the case of NbB 2 the phonon modes with peaks 
at 32, 81 and 106 meV couple to the electrons with almost equal strength, 
albeit much smaller than in MgB 2 , as can be seen from Fig. 2. We could 
have expected this because of the significant presence of the Nb d electrons 
at the Fermi energy. In Table III we have listed the Hopfield parameter 77, the 
electron-phonon coupling constant A, and the various averages of the phonon 
frequencies for NbB 2 and MgB 2 . The values listed in Table III for MgB 2 are 
in good agreement with the corresponding results of Choi et al. [5]. 

To examine the superconducting transition temperature, if any, of NbB 2 and 
MgB 2 we have used the calculated Eliashberg function a 2 F(u) to solve nu- 
merically the isotropic gap equation [7,16], and the results are shown in Fig. 
3 for a range of values of /i*. From Fig. 3 we find that for //* = 0.1 the T c 
for MgB 2 is equal to ~ 23 K, while for NbB 2 it is equal to ~ 4 K. Thus, our 
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Table 3 

The calculated Hopfield parameter rj, the average electron-phonon coupling constant 
A, the root mean square < u 2 > x l 2 and the logarithmically averaged uj\ n phonon 
frequencies for MgB 2 and NbB^. 



alloy 


7] {mRy / a.u?) 


A 


< u? >V2 (K) 


Uln (K) 


MgB 2 


167 


0.59 


835 


768 


NbB 2 


203 


0.43 


669 


494 



calculation shows that NbB 2 is superconducting with a possible T c of around 
3 K. It is worthwhile to point out that to obtain a T c close to that of 39 K 
for MgB 2 , as found experimentally [1], one has to solve the anisotropic gap 
equation [5] . The need to solve the anisotropic gap equation for MgB 2 rather 
than the isotropic gap equation arises due to the highly anisotropic electron- 
phonon coupling A(k, k') [3-5] over the Fermi surface. In the case of NbB 2 the 
electron-phonon coupling is neither as strong nor as anisotropic, and thus the 
results obtained with the isotropic gap equation are reliable. 

As indicated earlier, the experiments show superconductivity in hole-doped 
Nb x B 2 under pressure and Boron-enriched NbB 2 . Hole-doping and Boron en- 
riching both lead to a relative increase in the Boron population at the Fermi 
energy and, probably, enhances the peak around 81 meV leading to an increase 
in T c . Of course a more quantitative investigation is needed to pinpoint the 
exact nature of changes in NbB 2 which lead to superconductivity. 

In conclusion, we have studied the electron-phonon interaction in NbB 2 and 
MgB 2 , using full-potential, density- functional-based methods in P6/mmm 
crystal structure. We have described our results in terms of (i) electronic 
structure, (ii) phonon density of states, (iii) Eliashberg function, and (iv) 
the solutions of the isotropic Eliashberg gap equation, which clearly show 
significant differences in the electron-phonon interaction in NbB 2 and MgB 2 . 
We find that the average electron-phonon coupling constant is equal to 0.59 for 
MgB 2 and 0.43 for NbB 2 , leading to superconducting transition temperature 
of around 22 K for MgB 2 and 3 K for NbB 2 . 
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